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We develop a new model of the fluctuation dynamo in which the magnetic field is confined to 
thin flux ropes advected by a multi-scale model of turbulence. Magnetic dissipation occurs only via 
reconnection of the flux ropes. This model can be viewed as an implementation of the asymptotic 
limit Rm — * oo for a continuous magnetic field, where magnetic dissipation is strongly localized to 
small regions of strong field gradients. We investigate the kinetic energy release into heat, mediated 
by the dynamo action, both in our model and by solving the induction equation with the same flow. 
We flnd that a flux rope dynamo is an order of magnitude more efficient at converting mechanical 
energy into heat. The probability density of the magnetic energy release in reconnections has a 
power-law form with the slope —3, consistent with the Solar corona heating by nanoflares. 

PACS numbers: 95.30.Qd, 47.65.Md, 52.35.Vd, 96.60.Iv, 96.60.qe 



The dynamo action, i.e., the amplification of magnetic 
field by the motion of an electrically conducting fluid 
(plasma), is the most likely explanation for astrophysical 
magnetic fields. Evolution of magnetic field B embedded 
in a flow at a velocity u is governed by 

— = V X (u X B) + £B, (1) 

where C is an operator describing magnetic dissipation. 
In rarefied plasmas, such as the Solar corona, hot gas in 
spiral and elliptical galaxies, galactic and accretion disc 
halos, and laboratory plasmas, an important (if not dom- 
inant) mechanism for the dissipation of magnetic field 
is the reconnection of magnetic lines rather than mag- 
netic diffusion the latter modeled with C = r/V^ (if 
J] = const). Discussions of dynamos often refer to mag- 
netic reconnection, but attempts to include any features 
specific of magnetic reconnection to dynamo models are 
very rare On the other hand, theories of magnetic 
reconnection rarely, if ever, refer to the dynamo action 
as a mechanism maintaining magnetic fields. This paper 
attempts to bridge the gap between the two major ar- 
eas of magnetohydrodynamics by developing a dynamo 
model explicitly incorporating magnetic reconnections. 

The nature of the dissipation mechanism is important 
for the dynamo action. For example, dynamo action with 
hyperdiffusion, C = — r/iV* (and with a helical u) has 
larger growth rate and stronger steady-state magnetic 
fields than a similar dynamo based on normal diffusion 
0. This is not surprising as the hyperdiffusion operator, 
having the Fourier dependence of fc^, rather than fc^ of 
the normal diffusion, has weaker magnetic dissipation at 
larger scales. The release of magnetic energy in smaller 
regions (and larger current densities) in hyperdiffusive 
dynamos may also lead to a higher rate of conversion of 
kinetic energy to heat via magnetic energy. Magnetic 
hyperdiffusion also appears in the context of continuous 



models of self-organized criticality in application to the 
heating of the Solar corona The aim of such models 
is to reproduce the observed frequency distribution of 
various flare energy diagnostics. 

Magnetic reconnections may have an even more ex- 
treme form of the dissipation operator than the hyperdif- 
fusion: here magnetic fields dissipate only when in close 
contact with each other, so that the Fourier transform of 
C can be expected to be negligible at all scales exceed- 
ing a certain reconnection length do. It is then natural 
to expect that dynamos based on reconnections (as op- 
posed to those involving magnetic diffusion) will exhibit 
faster growth of magnetic field, more intermittent spa- 
tial distribution and stronger plasma heating. In this 
paper we consider dynamo action based on direct mod- 
eling of magnetic reconnection. For this purpose, we fol- 
low the evolution of individual closed magnetic loops in 
a model of turbulent flow (known to be a dynamo) and 
reconnect them directly whenever their segments come 
into sufficiently close contact, with appropriate magnetic 
field directions. As we show here, our model exhibits 
a power-law probability distribution of the magnetic en- 
ergy release similar to that observed in the Solar corona. 

Magnetic reconnection is usually modeled with the in- 
duction equation, £ = TyV^ (perhaps including the Hall 
current), and magnetic dissipation is enhanced due to 
the development of small-scale motions and magnetic 
fields. This approach may or may not apply to mag- 
netic fields concentrated into flux ropes, where magnetic 
energy losses are strongly reduced at large scales and, 
hence, more energy can be deposited at the smaller scale 
of order the tube radius, where reconnections occur. Our 
model explores this possibility. Furthermore, our model 
can be viewed as a numerical implementation of the limit 
i?ni — > oo for a continuous magnetic field, where magnetic 
dissipation is confined to strongly localized regions with 
exceptionally high magnetic field gradients. 




FIG. 1: The algorithm for inserting new trace particles in a 
stretched magnetic flux tube. If the distance between any two 
particles (shown with open circles) exceeds a length scale d, a 
new one is inserted between them, shown with a filled circle. 
Labels represent magnetic field strength. 
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FIG. 2: Reconnection occurs when the distance between two 
trace particles reduces to do (left); the connection of the par- 
ticles on a magnetic fiux tube changes after the reconnection, 
and the two closest particles are removed (right). 



We model the evolution of thin flux tubes, frozen into 
a flow, each with constant magnetic flux ijj. In this paper, 
we focus on the kinematic behavior, where the velocity 
field is independent of magnetic field. To ensure that 
V • B = 0, we require that our tubes always take the 
form of closed loops. Numerically, we discretize the loops 
into fluid particles and track their position and relative 
order (i.e., magnetic fleld direction) by introducing a flag 
denoted P, with P increasing along a given magnetic flux 
tube. Initially the particles are set a small distance apart, 
0.75c?, where d is an arbitrary (small) constant length 
scale. If, during the evolution of the loops, the distance 
between neighboring fluid particles on a loop becomes 
larger than d, we introduce a new particle between them, 
as illustrated in Fig. [TJ We use linear interpolation to 
place the new particle halfway between the old ones. The 
new separation between the particles is thus greater than 
0.5c? - this will be important when we consider removing 
particles. Thus, the spatial resolution of our model is d. 

Each particle is also assigned a flag B (Fig. [1]) for the 
strength of magnetic fleld at that point on the loop. As- 
suming magnetic flux conservation and incompressabil- 
ity, magnetic field strength in the fiux tube is propor- 
tional to its length. Magnetic fleld is initially constant at 
all particles, B = 1. When a new particle is introduced, 
magnetic field is doubled, as shown in Fig.[Tl at two out of 
three particles involved: this prescription emerged from 
our experimentation with various schemes and allows us 
to reproduce the evolution of magnetic fleld strength in 
a shear flow. Conversely, when the flow reduces the sep- 
aration of particles to less than G.5d, we remove a parti- 
cle. The value of the magnetic field strength fiag is also 
halved on the remaining particles in a manner consistent 
with the above algorithm. We have verified that this pre- 
scription reproduces accurately an exact solution of the 
induction equation for a simple shear fiow. 

If the separation between two particles, which are not 
neighbours, becomes less than a certain scale do, we re- 
connect their flux tubes by reassigning the flags P (Fig. (2) 
which identify the particles ahead and those behind of 
those involved in the reconnection. (To obtain meaning- 



ful numerical results, do has to be comparable to c?, e.g., 
c?o = l.5d.) Two particles are removed from the sys- 
tem after each reconnection event (and their magnetic 
energy is lost, presumably to heat). We also monitor the 
cross product of magnetic flelds close to the reconnection 
point. By ensuring that its magnitude is smaller than 
some tolerance e « 10~^ and that the magnetic flelds in 
the reconnecting loops are (almost) oppositely directed, 
we prevent parallel flux tubes with the same field di- 
rection from reconnecting. We monitor the amount of 
magnetic energy released in each reconnection event. To 
place the reconnection-based dynamo into a proper per- 
spective, we compare it with a dynamo obtained for the 
same velocity field, but by solving the induction equa- 
tion, i.e., Eq. IT]) with C = TyV^. In particular, we com- 
pare the rates of magnetic energy dissipation, which can 
be identified with the plasma heating rate. We assume 
that the part of the magnetic energy which drives plasma 
motion at a reconnection site (such as jets) is eventually 
dissipated into heat as well, so that we consider that the 
whole magnetic energy released is converted into heat. 
For the induction equation, the relevant quantity is 



d\nM/dt= / r/B-V^Bdl/ / dV 
Jv Uv 



(2) 



where M is the total magnetic energy. A similar quantity 
can be obtained for the reconnection-based dynamo by 
adding the contributions of all reconnection events to the 
magnetic energy release: 



7r = 



d\nM 
dt 



1 



SttMt ^ 

i—1 



BfSiLi 



(3) 



where r is a time interval during which Nt^ reconnections 
occur (we take r to be equal to ten time steps; indi- 
vidual reconnection events occur in a single time step), 
and Bi, Si and Li are the magnetic field strength, the 
cross-sectional area and length of the reconnected (and 
thus removed) fiux tube segment associated with a par- 
ticle number i. From our assumption of frozen flux. 
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BiSi = -0 = const, the total magnetic energy M is, 



^ Stt Stt ^ ' 

1=1 1=1 



(4) 



where A'tot is the total number of particles, and 

-1 



i=l 



'Nu 



y^ BiLi 



i=l 



(5) 



Any comparison of the solutions of the induction equa- 
tion with those from the reconnection model is not 
straightforward because of the difference in the control 
parameters of the two models: the magnetic Reynolds 
number = uolo/rj and the reconnection length do, 
respectively. A proxy for the magnetic Reynolds num- 
ber can be constructed from do as Rm = uqIq / (utdo) , 
where Mr is the characteristic reconnection speed. The 
reconnection-based dynamo is significantly more efficient 
than the hydromagnetic dynamo, in the sense that the 
growth rate of magnetic field in the former is significantly 
larger when R-^ k, Therefore, in order to achieve 

conservative conclusions, we compare dynamos with sim- 
ilar growth rates of magnetic field. Thus, i?ni > in 
the models compared. Magnetic field growth in a dynamo 
is obtained from the difference between the magnetic 
stretching and dissipation rates. In the reconnection- 
based dynamo, both are larger than those in a similar 
diffusion-based dynamo, but their difference is kept the 
same in the models which we compare below. 

We consider dynamos driven by two types of flow. 
Firstly, this is the Kinematic Simulation (KS) model of 
a turbulent flow 5], known to be a dynamo ^6(]. Here 
velocity at a position x and time t is 



u(x,t) 



N 

E 

n=l 



(A„ X k„ COS 0„ -I- B„ X k„ sin (/)„), (6) 



where 4'n = k„-x+w„i, N is the number of modes, k„ and 
'-^n — knUn are their wave vectors and frequencies. An 
advantage of using this flow is that the energy spectrum, 
E{kn) is controllable via appropriate choice of A„ and 
B„. We also note that V • u = 0. We adopt an energy 
spectrum which reduces to E{k) oc k~P for 1 ^ fc ^ k^, 
with k — 1 at the integral scale; p = 5/3 produces the 
Kolmogorov spectrum, and fc^r is the cut-off scale. We 
have adapted ([6]) to periodic boundary conditions. 
We also used the ABC flow of the form [7,] 

u — (cosy + sin z, sin a; -I- cos z, cos a; -I- sinj/) , (7) 

also known to support dynamo action, to demonstrate 
that our results are not sensitive to the form of the flow. 

The initial condition is a random ensemble of closed 
magnetic loops, and both the induction equation and the 
flux rope model are evolved with the same velocity field 
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FIG. 3: Magnetic energy release rates from the induction 
equation (dashed) and the flux rope model (solid) . The former 
has a mean value of 2.4 (here Rm = 1200) once the eigensolu- 
tion has developed. The latter (with Tim = 174) has a mean 
value of 23 (thick horizontal line). 



(apart from the overall normalization to provide compa- 
rable growth rates of magnetic field) . The initial condi- 
tion for the induction equation is obtained by Gaussian 
smoothing of the magnetic field in the ropes (this pro- 
cedure preserves V ■ B = 0). To evolve the induction 
equation, we use the Pencil Code [3] on a 256'^ mesh 
with 1000 < i?in < 1500 in a periodic box. The test 
particles in the fiux ropes are evolved using a 4'^ order 
Runge-Kutta scheme, with a time step of Zjv/(20ujv). 
The algorithm for inserting and removing points is ap- 
plied every time step, and the reconnection algorithm, 
every ten time steps. We choose c? to be 1/4 of the small- 
est length scale in the flow and set do/d = 1.5. 

Figure [3] shows the energy release rates in simulations 
where the growth rate of the magnetic field is cr = 0.16 in 
both simulations (with the unit time Iq/uq). The dashed 
line shows the energy release rate from a simulation of 
induction equation with R^^ = 1200, which has the mean 
energy release rate 71 « 2.4. The solid line shows the 
corresponding results from the fiux rope dynamo, with 
the mean value plotted as a dashed horizontal line. The 
mean value of the energy release rate from the reconnect- 
ing flux rope dynamo is 7r ~ 23, an order of magnitude 
larger. Also note strong fluctuations in the energy re- 
lease rate from the reconnection model, which are absent 
in the solutions of the induction equation. 

Dynamos with the ABC flow behave similarly. With 
i?m = 55, the induction equation gives an energy release 
rate of about 7i = 0.6. The corresponding flux rope 
dynamo with the same growth rate (0.02) has the energy 
release rate of 71 ~ 6.7, again ten times larger. 

Our approach is deliberately oversimplified with re- 
spect to the (incompletely understood) physics of mag- 
netic reconnection. Nevertheless, we can argue that our 
model is conservative with respect to the reconnection ef- 
ficiency. The reconnecting segments of magnetic lines in 
our model approach each other at a speed ~ uoRe""'^^^ 
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FIG. 4: Probability density for the normalized magnetic en- 
ergy release in individual reconnection events, AM/S^ms, 
from the time series of Fig. |3l for the flux rope dynamo (cir- 
cles) and the diffusive dynamo with the same magnetic fleld 
growth rate and velocity field form (squares). A power-law fit 
to the former and a Gaussian flt to the latter are shown solid 
and dashed, respectively. 



for the Kolmogorov spectrum, equal to velocity at the 
small scale do Ig with Iq the energy-range scale of the 
flow and do assumed to be close to the turbulent cut- 
off scale. If magnetic field is strong enough, the Alfven 
speed Va, which controls magnetic reconnection in more 
realistic models, is of order u{lo)- Then u^- <C Va and our 
model is likely to underestimate the efhciency of recon- 
nections. The Sweet-Parker reconnection proceeds at a 
speed of order VA-Rm^^^, whereas the Petschek reconnec- 
tion speed is comparable to VA/lni?m [H- For uq ~ Va 
and Rm — Re ^ 1, the reconnection rate in our model is 
larger than the former but much smaller than the latter. 

A remarkable feature of the energy release in the rope 
dynamo is that its probability distribution has a power 
law as shown in Fig. 21 /(x) oc where x = AM/B^^^^^ 
is the magnetic energy released in a reconnection event 
normalized to the mean magnetic energy, with the slope 
s ~ 3.3. Importantly, the same scaling, s ~ 3.0, 
emerges when we use the ABC flow instead of KS. A 
similar exponent arises in a reconnection model for the 
corona Q where, however, dynamo action is not included. 
Thus, weak 'flares' dominate the energy release in our 
reconnection-based system, as in the nanoflare model of 
coronal heating Interestingly more recent results 

with a nonlinear adaptation of the model retains this 
feature with s « —3.1 for the KS flow in the statistically 
steady state. We stress that the power-law behavior is 
not related to the self-similar nature of the velocity field: 
solution of the induction equation with the same velocity 
field, also shown in Fig. [U has an approximately Gaus- 
sian probability distribution. It is not as yet clear if the 
flux rope dynamo represents a physical example of self- 
organized criticality, but the system does possess some of 
the required properties. In particular, our reconnection 
model has a natural threshold in terms of the current 
density J > Bmin/'^Oi where iJmin = 1 is the minimum 
magnetic field, and, as we argue above, our model can be 
viewed as an extreme case of magnetic hyperdiffusivity. 



Furthermore, our simulations are kinematic (so, magnetic 
energy density is assumed to be small), whereas the So- 
lar corona is magnetically dominated. The importance 
of this distinction needs to be carefully investigated. 

To summarize, we have confirmed that the dynamo ac- 
tion is sensitive to the nature of magnetic dissipation and 
demonstrated that magnetic rcconnections (as opposed 
to magnetic diffusion) can significantly enhance the dy- 
namo action. We have explored the kinematic stage of 
the fiuctuation dynamo in a chaotic fiow that models 
hydrodynamic turbulence and in the ABC flow, with the 
only magnetic dissipation mechanism being the reconnec- 
tion of magnetic lines implemented in a direct manner. 
In our model, where magnetic dissipation is suppressed 
at all scales exceeding a certain scale do , the growth rate 
of magnetic field exceeds that of the fluctuation dynamo, 
based on magnetic diffusion, with the same velocity fleld. 
Even when the velocity fleld of the reconnection-based 
dynamo is reduced in magnitude as to achieve similar 
growth rates of magnetic energy density, the rate of con- 
version of magnetic energy into heat in the reconnection 
dynamo is a order of magnitude larger than in the corre- 
sponding diffusion-based dynamo. Thus, rcconnections 
more efficiently convert the kinetic energy of the plasma 
flow into heat, in our case with the mediation of the dy- 
namo action. This result, here obtained for a kinematic 
dynamo, can have serious implications for the heating of 
rarefled, hot plasmas where magnetic rcconnections dom- 
inate over magnetic diffusion (such as the corona of the 
Sun and star, galaxies and accretion discs). In contrast to 
the fluctuation dynamo based on magnetic diffusion, the 
probability distribution function of the energy released in 
the flux rope dynamo has a power law form not dissimilar 
to that observed for the Solar flares. 
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